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ABSTRACT 

We present the first complete 12 CO 7 = 3-2 map of M81, observed as part of the Nearby Galaxies Legacy Survey being carried out at 
the James Clerk Maxwell Telescope. We detect nine regions of significant CO emission located at different positions within the spiral 
arms, and confirm that the global CO emission in the galaxy is low. We combine these data with a new Ha map obtained using the 
Isaac Newton Telescope and archival H i, 24 //m, and FUV images to uncover a correlation between the molecular gas and star forming 
regions in M81. For the nine regions detected in CO 7 = 3 - 2, we combine our CO 7 = 3-2 data with existing CO 7=1-0 data to 
calculate line ratios. We find that the ratio 7 = (3 - 2)/(l - 0) is in agreement with the range of typical values found in the literature 
(0.2 - 0.8). Making reasonable assumptions, this allows us to constrain the hydrogen density to the range (10 3 - 10 4 ) cm -3 . We also 
estimated the amount of hydrogen produced in photo-dissociation regions near the locations where CO 7 = 3-2 was detected. 

Key words, galaxies: individual: M81 - Galaxies: ISM - Galaxies: star formation - ISM: molecules - ISM: H n regions - ISM: atoms 



1. Introduction 

The details of the interplay between the gas content of galaxies 
and the processes triggering star formation form one of the main 
open questions in astronomy (the various reviews and panel dis- 
cussion transcripts in Knapen et al. 2008 give a broad overview 
of the topic). Although it is known that star formation occurs 
mainly in the interior of giant molecular clouds (GMCs), the 
role played by the different phases of the gas is still poorly 
understood. For example, Kennicutt et al. (2007]) demonstrated 
that molecular gas traces star formation significantly better than 
atomic gas. However, the lack of a permanent dipole moment 
makes the direct detection of molecular hydrogen difficult. This 
forces us to use tracers such as CO to calculate the amount and 
distribution of molecular gas. The conversion factor from CO to 
H 2 (the so-called X C o-h 2 factor) has a value that must be con- 
sidered carefully and the relation of which with the morpholog- 
ical type, metalli city, and other parameters of galaxies is still 
disputed (e.g., |Maloney & Black|p88||Wilson|1995 



Tomisaka 2005 jBolatto et al.|2008l|Pineda et al.|2008t 



2009] T 



Wada & 



Zhu et al. 
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On the other hand, the idea that star formation is a one-way 
process from gas to stars is an exceedingly simplistic scenario. 
It has been argued (see, e.g., Allen et al.|1985|[T997[|Smith et al. 
2000] |Heiner et al.|2008a|b| ) that at least part of the atomic gas 
observed could be produced by photodissociation of molecular 
hydrogen. In regions with strong star formation, the young, hot 
stars produce vast amounts of UV ionizing photons that can dis- 
sociate an important part of the molecular gas in the host GMC. 
Therefore, the ability to trace the distribution of molecular gas 
with good spatial resolution is of vital importance if we wish to 
understand the behaviour of all the factors involved in the star 
formation process. 

This paper intends to shed further light on these issues by 
presenting the first complete 12 CO 7 = 3-2 map covering the 
inner part of the optical disk of M81, obtained as part of the 
Nearby Galaxies Legacy Survey (NGLS) on the James Clerk 
Maxw ell Telescope (JC MT ). Several authors (see, for exam 



2009 



and 



Iono et al. 



pie, [Wilson et al.| 

7 = 3-2, being more energetic and of higher critical density 



2009fr show that 12 CO 



(n c ~ 2.5 x 10 4 cm" J ) than lz CO 7 = 1 - (n c ~ 700cm" J ), 
is a better tracer of the warm, dense gas involved in star for- 
mation. Thus, this new CO map, in conjunction with new high- 
resolution Ha data from the Isaac Newton Telescope (INT) and 
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Table 1. Summary of M81 properties. 



Parameter 



Value* 



Object name 
R.A. (J2000.0) 
Dec. (J2000.0) 
B magnitude 
Classification 

Heliocentric velocity (kms -1 ) 
Distance (Mpc) 
D 2 5 (arcmin) 
12 + log(0/H) 
Scale (pc/arcsec) 



M81.NGC 3031, UGC 05318 

9:55:33.5 

+69:04:00 

7.89 

SA(s)ab z? 
34 + 4 
3.64+0.34 c 
26.9 x UA d 
8.41 e 
17.5 C 



Notes. 



Data from de Vaucouleurs et al. 



j except where in- 

dicated otherwise {b) Data from fButa et al J 420071 | c) |Freedman et al.| 
( | 1994} {d) Major and minor axis {e) Moustakas]( [2006t 



(1991 



public FUV, H i, and IR images, offers the possibility to study 
in detail the interplay between atomic and warm molecular gas 
and their relation with star formation and photodissociation pro- 
cesses. 

The JCMT is currently hosting a series of Legacy Surveys, 
including the NGLS that aims to observe 155 nearby galaxies 
(within 25 Mpc). The sample is selected to cover a wide spec- 
trum of morphologies and gas content and overlaps partially with 
the Spitzer Infrared Nearby Galaxies Survey (SINGS ; |Kennicutt| 
et al.||2003 ) sample. The project entails two different phases: a 
first set of observations to obtain 12 CO / = 3 - 2 observations 
for the whole sample has already been completed, while the sec- 
ond phase, involving SCUBA-2 continuum observations at 450 
and 850 /mi, will start imminently. Further information about the 
survey and its science goals can be found in Wil son etaL| ( |2009| ). 

M81 is a nea rby (3.63 Mpc, [Freedman et al.||1994| ) nearly 
face-on SA(s)ab (Buta et al. 2007| ) spiral galaxy. The spatial 
scale is 17.5 pc per arcsec, which allows us a detailed study of its 
GMCs ( which have typical size s between 30 and 50 pc; see, for 
example, |Liszt & Burt on|1981| ). Some basic properties of M81 
can be found in Table [l[ Although M81 has been extensively 
studied at multiple w avelengths, just one c omplete CO map has 
been released to date. Brouill eTet al . | (p~9 88) observed M81 in the 
/=l-0 transition, finding weak CO emission across the galaxy. 
Subsequent observations in the same transition ( Brouillet et al. 
199T||Sage & Westpfahl|l^|Sage|1993||Sakamoto et al.l2001| 



Heifer et al.|2003HKnapen et al.|2006HCasasola et al. 



2007) have 



focused on small, well-defined regions, and confirm the apparent 
lack of strong molecular emission in M81. 

We begin by introducing (in Section [2) the observations and 
data processing employed in this study. A description of the main 
results obtained, including a detailed comparison of the selected 
regions through different wavelengths, is given in Section|3]and a 
co mparison with the ph oto-dissociation regions (PDRs) method 
of |Heiner et al. | ( [20 08a ) in Section]?] We calculate line ratios and 
gas masses in Section [5] Conclusions are presented in Section [6] 



2. Observations and data reduction 

2.1. CO J =3-2 map 

The 12 CO 7 = 3-2 data presented here were obtained as part 
of the JCMT NGLS. The observations were taken in raster map 
scanning mode during three nights in February 2008 using the 
16 receptors heterodyne array receiver (H ARP) with the au to- 
correlation spectrometer imaging system (Smith et al. 2003 ) as 



backend ( Wilso n et al. 2009 ; Warr en et al. 2010 ) The spectrome- 
ter was configured to achieve a bandwidth of 1 GHz with a reso- 
lution of 0.488 MHz (0.43 km s" 1 for the studied 12 CO 7 = 3-2 
transition). The JCMT beam at this wavelength is 14.5 arcsec 
FWHM and the signal was sampled using 7.276 arcsec pixels to 
satisfy the Nyquist criterion. Given the 17.5 arcsec spatial scale, 
this resolution allows us to resolve features separated ~ 250 pc 
from each other. The observed region for M8 1 covers a rectan- 
gular area of 13.2x7.7 arcmin, including the brightest regions in 
Ha, IR, and FUV (see Figure [T]). Following the criteria used for 
the whole NGLS ( Wilson et al. 2009 ), the galaxy was sampled 
until the target noise parameters (T* A = 19 mK at a resolution of 
20 km s" 1 ) were reached. 



The data reduction process is expl ained in detail in |Wilson 
|et al.| ( [2009] ) and in Appendix A of |Warren et aL] p010| ), so 
only a brief summary of the given steps is presented here. After 
flagging bad baselines and spikes, the individual scans were 
combined using a sine (nx) sine (knx) weighting function. The fi- 
nal cube was trimmed to remove the edges of the scans (ap- 
proximately 60 arcsec) where coverage is incomplete. To re- 
move baselines, line regions were masked out (applying differ- 
ent methods for narrow and broad lines and combining the two 
resulting masks) and third order polynomials were fitted to the 
continuum. These baselines were then subtracted from the orig- 
inal spectra. 

To compile the integrated CO map, we used a two-step 
approach. In the firs t step, we used the clumpfind algorithm 
( Williams et al. 1994 ) to select regions with emission above three 
times the RMS noise level (3<x), based on the contour levels de- 
fined by clumpfind in both the physical and spectral dimensions 
and on the number of channels in the data cube and the noise in 
each channel. This RMS value ensures, based on our experience, 
that all the potential features present in the CO map are found. 
We then applie d these masks to pro duce moment maps from the 
original cube ( Warr en et al.||2010| ). The data values were con- 
verted to main beam temperature by dividing by t]mb =0.6. 

In the second step, we imposed a different and more restric- 
tive criterion to each region in the moment map resulting from 
the first step, aiming to reach a detailed confirmation of every 
individual feature in the clumpfind CO map. We considered a 
region to be a significant detection only if the peak value of the 
integrated spectrum of the region was above 2.5<x, and if its ve- 
locity agreed with that expected from an H i velocity field. For 
the latter, we compared the CO spectra at each position with an 
Hi profile from THINGS ( [Walter et al.|2008] ) Hi data cube (see 
Section |2.3| for further details about these data). For each CO 
position, we compared the peak of CO emission with that at the 
equivalent position in the H i map, after smoothing both spectra 
to a channel velocity of 10 kms -1 . Figure [5] shows the CO and 
H i spectra for each detected feature. The agreement between the 
spectra, in terms of peak velocity and line width, is good in al- 
most every case. The average separation between the CO and 
Hi peaks is 4.4 ± 3.3 kms -1 . One spurious feature was rejected 
because its line velocity was in complete disagreement with that 
of neighbouring CO regions and with the expected velocity as 
derived from published Hi data (see below). 

As a result of this second step, a total of four features were 
manually rejected as their line emission peak was confirmed to 
be significantly below 2.5<x, or their line velocities or widths 
were differen t from the expected values from CO J =1-0 
and H i data (Hei fer et al.|[2003l |Knapen et aL][2006l |de Blok 
|et al.||2008| ). After these verifications, we confirm the detection 
of nine peaks of 2.5cr or more in M81. The final total intensity 
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(moment 0) map for M81 is displayed, along with the other maps 
used in this paper, in Figure [T] 

To quantify the uncertainty in our profiles, we estimated the 
noise level (by masking the CO line and calculating the standard 
deviation of the remaining spectrum) and compared it with the 
value of the peak. The results, shown in Table |2j confirm that all 
our detections are around 3<x or better. The noise levels calcu- 
lated are in good agreement with the global values measured in 
other galaxies of the NG LS ( [Wilson et aL]|2009[ |Warren et aT 
50T0l |Bendo et al.1[20TQ| ). Table |2j also includes the integrated 
flux ( J r MB dv) and the FWHM for each spectral line derived 
from a Gaussian fit. 



2.2. R and Ha images 

We present new R and Ha images acquired using the Isaac 
Newton Telescope in the course of a programme that aims to ob- 
tain Ha continuum- subtracted images for all the JCMT NGLS 
galaxies. The Wide Field Camera (WFC) was used to obtain 
three Ha exposures of 480 s each and two 7?-band exposures of 
120 s each. The WFC consists of four 2048x4100 pixel CCDs 
with 1 arcmin separation between the chips and a pixel scale of 
0.33 arcsec per pixel. The final mosaic covers an area of approx- 
imately 33x21 arcmin in right ascension and declination respec- 
tively. Although maps of similar resolution have been released in 
the past (for example, as part of the ancill ary data for the Spitze r 
infrared nearby galaxies survey, SINGS, Kenni cutt et al.|2003] ), 
our new map covers a larger area and is deeper than other maps 
of M81 published to date. The spatial resolution of the images 
is ~ 1.2 arcsec FWHM, which will allow us to detect features of 
the order of the typical size of GMCs. 

The reduction was carried out using IRAlQ following the 



prescriptions given in Knapen et al.| ( 2004] ) adapted for the WFC. 



After the bias and flat field correction, the images were normal- 
ized by dividing by the exposure time and then aligned and com- 
bined using a median rejection algorithm in order to remove cos- 
mic rays from the final image. The background level was calcu- 
lated by taking an area with little emission, rejecting those pixels 
with values 3<x above or below the mean, and iterating until con- 
vergence was achieved. This value was then subtracted from the 
image. Some defects caused by scratches on the CCDs were cor- 
rected, replacing the affected pixels with a linear interpolation of 
the nearest good pixels. 

At this point, we aligned the final Ha and 7?-band 
background- subtracted images. To subtract the continuum emis- 
sion from the Ha map, we plotted the intensity, in counts per 
second, of each pixel in the Ha image versus the intensity of the 
same pixel in the 7?-band image. Three percent of the brightest 
pixels were manually rejected to avoid the contribution of fore- 
ground stars. The resulting plot should display a straight line in 
the absence of Ha emission, and the slope of the line fitted to 
the upper limit of the distribution provides the scaling factor be- 
tween both images (see Figure 1 in Knape rTet al.|2004 ). Using 
this scale factor, the 7?-band continuum image was scaled and 
subtracted from the Ha image. The result is shown in Figure [T] 

2.3. Ancillary data 

We u sed the VL A Hi data cube obtained by B. Hine and A. 
Rots ( Hi ne|1984|). This fi nal column density map is the same as 
that used in | Allen et aL| ( |1997| ) and the details concerning the 
calculation of the moment map can be found there. The spatial 



resolution of the image is 9 arcsec. A new deeper Hi moment 
map of M81, based on archival VLA data, was released by the 
THINGS team ( [Walter et al.|2QQ8j ). The resolution of these new 
data is slightly better: 7.58x7.45 arcsec. Unfortunately, this map 
is not continuum subtracted and cannot be used for our compar- 
ison with other wavelengths. However, as described in Section 
2.1 we use the THINGS data to compare the position and width 
of the spectral lines found in the CO map. The presence of con- 
tinuum emission is not a problem for this use. 

The FUV image, obtained from the GALEX public 
databas^] was acquired as part of the G ALEX nearby galaxy 
survey (NGS; |Gil de Paz et al. 2004). The spectral range 
observed is 1344 - 1786 A w ith an effective wavelength of 
1528 A ( [Morrissey et al.||2007| ). The FOV is 1.28 degrees with 
a FWHM of 4 arcsec. The data can be converted into UV flux 
(erg cm -2 s" 1 A" 1 arcsec" 2 ) by ap plying a conversion factor of 
6.22 • 10" 16 ( [Morrissey et al.|2QQ7| ). 

We used the 24 micron image from SINGS to trace dust 
emission associated with recent star formation activity (Calzetti 



|et~aLl [20051 [20071 |Prescott et aLl[2007] ). The data were taken 
using the Multiband Imaging Photometer for Spitzer (MIPS; 
Rieke et al.| | 2004] ) in scan map mode. The final image is 



32.5 x 55 arcmin. The data have a FWHM of 6 arcsec (Spitzer 
Observers Manual), and the flux calibration accuracy is 4% 
( [Engelbracht et al.|2007| . The survey websit^] gives further de- 
tails on the data processing. 

For the Ha, 7?-band, FUV and 24yum data, the astrometry 
of the images was checked by comparing the positions of fore- 
ground stars with the Guide Star Catalogue ( [Lasker et al.|1990 ). 
It was found that the astrometry of all the images was precise to 
within 1 arcsec, which is sufficient for the purposes of this work. 



2.4. CO J =1-0 and J = 2-1 maps 

Although CO J = 1 -0 observations of t he inner regions 

( [Sakamoto et al.||2001HCasasola et al.||2007| ) and the western 
arm of M81 ( [Knapen et al.|2006| ) have been performed, the only 
comprehensive C O J = 1-0 map a vailable for M81 is the 
one published by Brouill eTet al. | ( T988] ) (hereafter the Brouillet 
map) using 12 meter National Radio Astronomy Observatory 
(NRAO) data. The half power beam width (HPBW) of the map 
is 60 arcsec. This map was later improved with 23 arcsec HPBW 
IRA]V03Om observations in some regions ( [Brouillet et al.|l99l| ). 
However, there is no spatial correspondence between their detec- 
tions and our CO J = 3 - 2 emission, thus the IRAM map cannot 
be used for further comparisons with our data. 

ISakamoto et aL| ( |200 1) , in the CO / = 1 -0 line, and |Casasola[ 
et al.|fl2007|), in both the J = 1 - and J = 2 - 1 lines, probed 



a sma ll area in the nu clear region of M81 (40 arcsec, Cas asola 
|et al.| and 20 arcsec, Sakamot o et al.) . We convolved our data 
cube to the resolution of their data and found no CO 7 = 3-2 
emission in those regions. This is no surprise as, at our current 
RMS of 19 mK at 20kms _1 resolution, their lines (of typical 
temperature ~ 50 mK, which with a typical CO J = 1 - to 3-2 
ratio would correspond to 7mb(3 - 2) ~ 25 mK) are too weak to 
allow CO J = 3 - 2 detections at a level of 2.5<x. 

|Knapen et al. ([2006 ) used the 45 m radio telescope at the 
Nobeyama Radio Observatory (NRO) to observe a small region 
in the western spiral arm of M8 1 . Unfortunately, this region is 
outside our mapped area. 



http ://iraf .noao . edu/ 



2 http : //www. galex . caltech. edu/index. html 

3 http://sings.stsci.edu/ 

4 Institut de Radioastronomie Millimetrique 
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Table 2. CO spectral analysis. 



Position 


RA a 


Dec fl 


Peak velocity 


Peak value (T MB ) 


cr b 


Peak-cr ratio 


FWHM C 


J r MB dv c 








kins" 1 


K 


K 




kms -1 


Kkms" 1 


1 


9:55:19.0 


+69:08:47 


139.5 


0.15 


0.041 


3.7 


15 


2.33 


2 


9:55:15.5 


+69:08:55 


156.4 


0.11 


0.039 


2.8 


6.5 


0.76 


3 


9:55:33.0 


+69:07:28 


88.7 


0.16 


0.047 


3.4 


8.5 


1.46 


4 


9:55:08.5 


+69:07:52 


164.9 


0.13 


0.043 


3.0 


17 


2.25 


5 


9:56:25.0 


+68:59:33 


-227.0 


0.11 


0.038 


3.0 


13 


1.55 


6 


9:56:04.0 


+68:59:01 


-238.0 


0.10 


0.035 


2.8 


28 


2.81 


7 


9:55:48.0 


+68:59:20 


-227.0 


0.14 


0.044 


3.1 


19 


2.73 


8 


9:55:37.5 


+68:59:04 


-190.6 


0.12 


0.042 


2.9 


8.2 


1.09 


9 


9:56:03.0 


+68:59:15 


-238.0 


0.12 


0.046 


2.7 


12 


1.56 



Notes. (a) Right ascension and declination of the peak pixel (J2000) (b) Calculated after masking +50 kms 1 around the line. Average cr is 0.042 
at 10 kms -1 resolution. (c) Obtained using Gaussian fits. 



2.5. Convolution kernels 

To compare images at different wavelengths it is necessary to 
smooth all the data to match the image with the lowest resolu- 
tion, which in this case is the CO J = 3 - 2 map (14.5 arcsec). 
For the Ha, R-band, and FUV images, we convolved the original 
data with gaussians of shape 

iff(r) oc e^/^co-^ 2 ), 

where cr is the original resolution of the image and cr co is the 
resolution of the CO map. The convolution was performed using 
the task gausmooth in Starlink/KAPPA. 

The MIPS 24 micron PSF cannot be approximated as a 
Gaussian function because of its strong Airy rings, so we can- 
not apply the above equation to smooth the data. Instead, we 
smoothed the data using a co nvolution kernel creat ed by follow- 
ing the methods outlined by |Gordon et al. ( 2008| ) and |Bendo| 
et al.| ( |2010|) using an empir ical 24 micron PSF originally de- 
rived by Young et al. ( 2009 ). The convolution was performe d 
using the task convolve in Starlink/KAPPA ( [Currie et aTf2 008 ). 

3. Detailed analysis of CO emission regions 

As presented in Figure[T] the emission of the star formation trac- 
ers (FUV, Ha, and 24 /mi) follows a common pattern with strong 
centrally peaked emission and well-traced spiral arms. This is 
also reflected in the H i map, except for the lack of nuclear emis- 
sion. All the detected CO J = 3 - 2 features are located in two 
small regions within the spiral arms. The comparison between 
Ha and FUV shows a good agreement. In almost every case, an 
individual region with enhanced Ha emission has a counterpart 
in the FUV map. The opposite is not always true. This result 
confirms those presented by Alle n et aL] (p"997). In this section, 
we compare in detail the emission of the different tracers on the 
smallest resolved spatial scales at the location of the CO emis- 
sion. 



Table 3. Statistics of the CO positions. 



3.1. CO regions 

Our CO emission is located in two regions of the galaxy: one in 
the northern and another in the southern arm. We call these re- 
gions A and B, respectively (see Figure [T]). Figures [3] and [4] show 
the detected CO features overlaid on the various other maps for 
each region. Although positions 1 and 2 (region A) are located 
within the same overall complex, we treat them as two different 
peaks; this is evident from the CO contour levels (Figure [3]). The 
same holds for positions 6 and 9 in region B (Figure [4]). 



Position 


FUV 


Ha 


24 //m 


Hi 


1 




C 




C 


2 


C 


C 


C 


C 


3 








C 


4 


C 




P 




5 










6 










7 










8 
9 


C 




c 




P(%) 








11 + 10 





C(%) 


33 + 16 


22+14 


33 + 16 


33 + 16 


RP (%) 


3 + 2 





1 + 1 


1 + 1 


RC (%) 


11 + 3 


10 + 3 


4 + 2 


15 + 4 



Notes. For each position and wavelength, the table indicates whether 
a peak of emission was found at the exact location (P) or if a peak is 
closer than 150pc (C). The two lines below indicate the percentage of 
detections in each case. The last two lines show the same statistics for 
the peak (RP) and close (RC) random positions. Note that the close 
percentages also includes the peak cases. Errors have been calculated 
using Poisson statistics. 



A detailed comparison with data at FUV, Ha, IR, and 
Hi wavelengths was carried out for each CO feature, with the 
results shown in Table [3] Although the typical size of a GMC is 
in the range of 30 - 50 pc, the accuracy in the determination of 
position and distances between features is limited by the FWHM 
of the CO map (14.5 arcsec, or -250 pc in M81). Taking these 
values into account, we consider two kinds of correlation: the 
case where a peak of the studied wavelength is found at 50 pc 
or less from the peak in the CO feature (marked as P, peak, in 
Table[3]), and the case when a peak is found in an area of 150 pc 
in radius around the peak of CO emission (marked as C, close). 
This task was performed for each wavelength using the images 
convolved to the CO resolution. 

A peak of emission close to the CO feature is found in all 
bands for 33% of the cases, except in Ha, where the probability 
is slightly lower (22%). On the other hand, the probability of 
finding an emission peak at the exact location of the CO peak 
is much lower (just one coincidence among all the positions and 
wavelengths), which indicates that the molecular gas does not 
occupy the same space as the atomic hydrogen or the ionized 
gas, on the scales considered. 
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3.2. Randomly selected regions 

The results presented above seem to indicate that the CO emis- 
sion is correlated with both Hi and star formation tracers. To 
place these results in perspective, we used a random number 
generator to obtain 100 completely random positions within the 
footprint of regions A and B. 

The last two rows of Table [3] present the results of the com- 
parison between these random regions and the emission at the 
available wavelengths. We find that the probability of a position 
lying exactly on a peak is quite low (3% for FUV, lower for the 
remaining wavelengths). When a peak is located near the central 
position (closer than 150 pc), the numbers do not change signifi- 
cantly: the probabilities increase to 15% only for Hi and always 
remain below the percentages obtained for the CO regions. 

These values, compared to the significantly higher ones 
found in the previous section for CO, indicate that a relation be- 
tween the CO J = 3-2 emission and the tracers of star formation 
and atomic gas indeed exists. The correlation is, however, not as 
strong as in the case of the star formation tracers and H i, where 
a corre spondence has bee n found for about 90% of the peak po- 
sitions ( [Alien et al.fT997) ). 



4. Analysis of candidate PDRs 

In principle, it is possible to estimate the total hydrogen vol- 
ume density of gas near a PDR from the measured UV flux and 
Hi column density, by making certain assumptions about the 
physical properties of the PDR. Comparing this estimate with 
measurements of the molecular gas column density, obtained 
from CO observations, one can then obtain information about 
the physics and even internal structure of the PDR, and constrain 
the CO to H2 conversion factor. 

Using the approach outlined in Hein er et aL| ( |2008a| ), we in- 
vestigate the presence of candidate PDRs near the detected CO 
J = 3 - 2 emission. While Hein er et~aL| selected the brightest 
FUV regions in the galaxy, we used our CO detections as ref- 
erence points to locate the nearby peaks of UV emission and 
Hi patches. 



The method applied was explained in depth in Heine r et al 
( 2008a ). However, for the sake of clarity, we include here a short 
explanation of the main steps. Allen] ( [5004 his equation 6) de- 
rived the following expression for the Hi column density pro- 
duced by photodissociation near UV sources: 



7.8 x 10 20 t 



1 + 



106 x Gn 



So 



-1/2 



This relation relies heavily on the dust-to-gas ratio (6/60, with 
£0 the local solar neighborhood value), which is derived from 
metallicity measurements, where Go is an expression of th e in- 
cident UV flux at the H 1 patch in units of the Habing flux ( |van| 
|Dishoeck & Black|[T988l |Menp004) and depends on the dis- 
tance to the galaxy and the separation between the UV source 
and the Hi column, pm. This expression can be inverted to cal- 
culate the total density of gas in the GMC, n(cm~ 3 ) = n^ + ln^ . 
At the core of the GMC , we can then expect that n = 2nu 2 . 

The parameters used here are the same as in Heiner et al.| 
(200 8a), except for the extinction [A FUV = 1.9xE(B- V) = 0.63 
([Gil de Paz et al.||2007] ) with E(B - V) = 0.080 (Schlegel| 
et al.|1998j >] and the calibration of the S /So relation [log(S/So) 



Table |4] shows the values measured for each of our CO de- 
tections and the total density of gas calculated for each region. 
Owing to the lower UV extinction value used in this work, the 
total density results are approxim ately half of those th at would 
be obtained with the values used in Heiner et al. (2008a). The un- 
certainties we re calculated as explained in Appendix A of|Heiner 



et al 



(2008a). The values obtained are in the range 1-178 cm 



-0.04 5 x/? ga i + 0.31, now using an updated value of solar metal- 
licity ( Allende Pri eto et aTf2 001 )]. We used the distance referred 
to in TableHJ 



for the total hydrogen volume density. The highest value occurs 
near CO position 2, which was already noted in Table [3] to have 
FUV, Ha, 24//m, and Hi emission within 150 pc, making it the 
candidate PDR most likely to be directly comparable to the CO 
emission at the same location. 

We speculate that the gas clouds traced by the PDRs are fun- 
damentally different from those detected in CO J = 3 - 2. The 
latter are indicative of a dense, hot environment, concentrated in 
a relatively small volume of space, where the PDR method is 
sensitive to lower density gas. A molecular cloud like this does 
not have to contain an OB association. For example, CO posi- 
tion 8 corresponds to a known H 11 region that is too faint to be 
detected in the GALEX image that we used. It might also corre- 
spond to a very dense core of a larger association of GMCs. An 
important additional clue to Hi being produced in PDRs, in turn 
produced by OB star clusters, would be the presence of bubbles 
of H 1 , where the atomic gas has been cleared by recent star for- 
mation. Comparing the locations of our CO detections with the 
Hi holes presented in Bagetak os et al. ( 2010 ^, we find that the 
majority of our detections occur on the edge of an Hi hole (the 
smallest Hi hole detected being 90arcsec across). Exceptions 
are positions 4, 7, and 8 where no H 1 holes were detected in the 
immediate vicinity. 

The gas densities derived with the PDR method are most 
likely to yield lower limits -we take the maximum H 1 column 
density, where the higher-density, more-compact PDRs might 
produce a lower H 1 column density because the column is devel- 
oped over a shorter distance. Those would be harder to discern 
in the H 1 map. We may also underestimate the incident UV flux, 
although not by much more than the GALEX detection limit in 
this case, so this effect would be small, but may be amplified if 
the separation between the UV source and the H 1 patch is small. 

Only near position 2 do we find, perhaps because of a 
favourable morphology of the gas, an HI patch sufficiently close 
to the UV source to trace a higher-density GMC. This position 
was not included in the candidate PDRs selected by Heiner et al. 
( 2008a ), who did not aim to study a complete set of FUV regions 
in M81. In contrast, we do not detect any of the GMCs traced 
with the PDR method in their paper, potentially because the low 
CO emission in those regions and M81 overall is insufficient to 
be traced in our map. 

5. Line ratios 

We calculated the CO (3 - 2)/(l - 0) line ratios for each of our 9 
CO / = 3 - 2 detections using the CO 7=1-0 measurements 
from |Brouillet et aL| ( |1988| ) derived from NRAO observations. 
To compare both data sets, we convolved our CO 7 = 3-2 
cube to match the resolution of the Brouillet data (60 arcsec; see 
Section |24| ). A Gaussian fit was then made to each line. In Table 
[5] we show the values and ratios for the line fluxes. At 60 arcsec 
resolution, positions 2 and 9 cannot be separated from 1 and 6, 
respectively, thus we give the average values. The Brouillet map 
does not cover position 8. 

The CO (3 - 2)/(l - 0) ratios for the integrated flux of each 
peak lie in the range (0.2 - 0.8) with a mean value of 0.43 ± 0. 1 8, 
which is consistent with previous results from the JCMT NGLS 
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Table 4. PDR data 



CO position 


FUV coordinate 


^gal 


TTV flnv 

u v nux 


Phi 


^Hl 


O/Oo 




// 


uncertainty in n 




K./\. 


Dec. 


kpc 


1A _15 _2 -1 * -1 

10 ergs cm s A 


pc 


10 cm 






-3 

cm 


% 


1 b,c 
1 


09:55:16.75 


+69:08:55.62 


J.O 


A 1 1 
4. 1 1 


4o>> 


Z.jU 


1 1 

1.1 jy 


U. 14 


n a 
U.4 


A A 
44 


2 h 


09:55:16.75 


+69:08:55.62 


5.6 


4.11 


26 


2.28 


1.139 


47.93 


178 


146 


3 


09:55:35.02 


+69:07:34.69 


4.7 


6.99 


81 


2.93 


1.251 


8.17 


7 


69 


4 


09:55:10.25 


+69:07:46.30 


4.7 


0.70 


182 


2.63 


1.260 


0.16 


0.2 


53 


5 


09:56:24.69 


+68:59:17.12 


8.0 


1.59 


251 


1.55 


0.887 


0.19 


5 


31 


6/9 d 


09:56:03.79 


+68:59:05.59 


6.0 


4.73 


278 


4.31 


1.100 


0.47 


0.1 


68 


1 


09:55:47.43 


+68:59:23.33 


5.1 


5.46 


95 


4.26 


1.205 


4.61 


0.6 


81 


8 e 


09:55:37.50 


+68:59:04.00 


6.0 


0.15 


170 


0.92 


1.096 


0.04 


1.8 


32 



Notes. See Section|4]for an explanation of the parameters in this table (a) Hi background subtracted was 10 20 cm -2 ib) CO positions 1 and 2 have 
one FUV source in the middle (c) A weak FUV source is present near position 1 but its flux could not be accurately measured (J) One measurement 
was carried out for CO positions 6 and 9 (e) Upper limit for a source at the GALEX detection limit (F FUV = 0.15 x 10~ 15 ergs cm -2 s _1 A -1 ). 



Table 5. Line ratios. 



CO J = 3 - 2 flux CO / = 1 - flux (3-2)/(l-0) 





Kkms" 1 


Kkms- 1 




1 


0.33 


0.8 


0.41 


2 


0.24 


0.8 


0.31 


3 


0.41 


0.5 


0.82 


4 


0.22 


0.8 


0.28 


5 


0.15 


0.7 


0.21 


6 


0.22 


0.5 


0.44 


7 


0.26 


0.5 


0.52 


8 


0.15 






9 


0.31 


0.7 


0.44 


Mean 






0.43 



( Wilson et al. 2009] |Warren et al.|2010HBendo et al.|201Q|[Irwin| 
et al.|2010| B. Tan et al. in prep.). 

Using th e JCMT data of N GC 4631 and a single compo- 
nent model, Irwi n et aL] ([2010) calculate the relation between 
the / = (3 - 2)/(l - 0) ratio and the hydrogen density for a 
number of temperatures. Assuming temperatures from 10 K to 
50 K and a CO (3 - 2)/(l - 0) ratio in the range 0.2 to 0.8, this 
model constrains the hydrogen density, n(H 2 ), in M81 to a range 
between 10 3 and 10 4 cm -3 . 

Israel (2009) proposed the use of a two component model as 
the best approach to studying the molecular gas dynamics. In this 
model, the warm (T k i n ~ 150 K) gas represents just a small part 
of the overall molecular gas content, which is mostly (~ 80%) 
cold (7ki n ~ 10 K). With this additional constraint, the hydrogen 
density for M81 is most likely to be near the lower end of the 
range indicated above. 



6. Concluding remarks 

We have presented a new CO 7 = 3-2 map of M8 1 obtained 
as part of the JCMT NGLS, as well as a new Ha image from the 
INT, and combined these new data with public H i , IR, and FUV 
data. The study of these data leads us to the following conclu- 
sions: 

(i) We confirm an overall low level of CO emission in M8 1 . The 
emission that we do detect is limited to a few small regions 
within D25/2, all of which are located in the spiral arm re- 
gion. These results agree qualitatively with those previously 
published for the / = 1-0 transition of CO. Moreover, the 
regions where the CO is detected also seem to be the areas 
where the star formation is more conspicuous (see Figure [T]). 



(ii) We compare the spectral profiles of our CO / = 3 - 2 de- 
tections with those of H 1 data at the same positions. There is 
good agreement between the peak velocities and line widths 
in both cases, which confirms our CO detections and indi- 
cates that the kinematical conditions of the molecular and 
atomic gas are similar. 

(iii) We confirm the previous results by Alle n et al.| ([T997) that, 
locally, the Ha emission is almost without exception well 
traced by the FUV emission. The opposite is not always true, 
and some features in the FUV map do not have any corre- 
sponding detection in Ha. The same holds when comparing 
Ha with the 24 jim emission. 

(iv) We found statistical evidence supporting the relation be- 
tween molecular gas and star formation in M81, although 
the correlation between the star formation tracers is much 
stronger. 

(v) We investigated candidate PDRs close to the CO 7 = 3-2 
detections and find total hydrogen volume densities in the 
range 1-178 cm" 3 , in a greement with other candidate-PDRs 
studied in |Heiner et al. ( 2008a| ). The morphology at the loca- 
tion of the highest density is indicative of a relatively small, 
compact PDR, as opposed to the other locations where such 
a morphology could not be seen. 

(vi) We calculated / = (3 - 2)/(l - 0) line ratios for each of our 
detections, obtaining values in the range (0.2 - 0.8) with a 
mean value of 0.43+0.18. This is consistent with previous 
work for other galaxies in the NGLS. Using a single compo- 
nent model we constrained the hydrogen density to the range 
(10 3 - 10 4 ) cm -3 , most likely near the lower value of this in- 
terval. 
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Fig. 2. CO J = 3 - 2 spectra (blue solid) obtained from the peak pixel of each position where we detected CO, compared to the 
H i spectra (red dashed) at those positions. In each case the spectra have been smoothed to Av = 10 km/s and the H i profile has been 
scaled to fit the CO peak (pointed out with an arrow at each position). No spatial smoothing was applied. The green dashed lines 
mark the cr ~ ±0.04 K average level (see Table [2]). The positions labeled are the same as in Figures [3] and [4] 
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